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ABSTRACT 

We use high-resolution simulations of large-scale structure formation to analyze the effects of inter- 
acting dark matter and dark energy on the evolution of the halo mass function. Using a chi-square 
likelihood analysis, we find significant differences in the mass function between models of coupled 
dark matter-dark energy and standard concordance cosmology (ACDM ) out to redshift z — 1.5. We 
also find a preliminary indication that the Dark Energy Survey should be able to distinguish these 
models from ACDM within its mass and redshift contraints. While we can distinguish the effects 
of these models from ACDM cosmologies with different fundamental parameters, DES will require 
independent measurements of a$ to confirm these effects. 

Subject headings: cosmology:theory, dark matter, dark energy, structure formation, methods: N-body 
simulations 



1. INTRODUCTION 

Multiple independent lines of evidence, including ob- 
servations of the large-scale matter distribution (eg. 
iPercival et al. l2001l). cosmic microw ave background fluc- 
tuations (eg. Hinshaw et al. 2008), and type la super- 



novae (jPerlmutter et all 11999k iRiess et all 11998ft , sug- 
gest that our universe is dominated by two compo- 
nents: dark matter, which is probably a form of non- 
baryonic matter, and dark energy, which is a name for 
the presently unknown cause of the observed acceleration 
of expansion. However, we still lack an understanding 
of any possible i nteractions betw een these two princi- 
al constituents jBean et aill2005[ ). In a previous paper 
Sutter fe R ickcr 20Q8J; hereafter SR08), we examined the 
role that interacting dark matter and dark energy would 
play in the development of one-dimensional Zel'dovich 
pancakes, an important idealized case useful for under- 
standing structure formation. Here we extend that pre- 
liminary work to a more realistic three-dimensional sim- 
ulation of the growth of dark matter halos in an attempt 
to find ways to distinguish these models from standard 
cosmology. 

In this paper, we study the effects of a Yukawa in- 
teraction between a single family of nonrelativistic dark 
matter (DM) particles and a scalar field that is respon- 
sible for the dark ene rgy (DE). We fo llow closely the for- 
malism developed bv iFarrar fc Peebles! ((2004). Such an 
interaction is initial ly attractive bec ause it is motivated 
by particle physics (| Amendolal l2000t ) and might provide 
a way to explain the apparen t emptiness of the voids, as 
demon strated numerically bv lNusser. Gubser. fc Peebles! 
(|2005h . 

There has been considerable interest recently in study- 
ing the effects of these int eractions on structure, both 
using an analytic approach (Mai nini fc Bonomettdl200a l 
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and using direc t simulations (eg. iMaccio et al.l l2004t 
iManera fc Mota 2006; and others). However, the cur- 
rent numerical studies suffer from poor resolution, and 
these results cannot reliably be compared to simulations 
of standard cosmological structure formation. In this pa- 
per, we use high spatial resolution and careful analysis 
to accurately capture many dark matter halos for use in 
comparison. 

The following is a brief summary of the equations we 
solve and our numerical techniques. In Section [2] we dis- 
cuss modifications to the halo mass function. We use 
these mass functions to distinguish interacting DM-DE 
from standard concordance cosmology using a x 2 like- 
lihood test. Additionally, we discuss the feasib ility of 
using the Dark Energy Survey (|Annis et al.ll2005h to de- 
tect this coupling within its mass and redshift contraints. 
Finally, we determine the extent to which two specific 
models of interacting DM-DE can be differentiated from 
each other. 

1.1. Analytical Methods 

Compared to simulations of the full non-linear theory, 
we found in SR08 that the perturbation theory presented 
by Farrar and Peebles is very accurate in determining 
the evolution of structure, and hence we will maintain 
the perturbative approach and assume fluctuations in the 
scalar field are small. Under perturbation theory, the ho- 
mogenous part of the dark energy scalar field, 4>b, evolves 
as 

a , o«i , dV , o n "ioH'o 1 . 



8ttG (j) b 



0. 



(1) 



where fl m o is the dark matter particle fraction of the crit- 
ical density and Hq is the Hubble constant. A subscript 
of denotes the present-day value. The dark matter par- 
ticle equation of motion is 
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TABLE 1 
Simulation potential parameter value choices. 

Here $ is the normal comoving gravitational potential, 
a is the scale factor, v is the comoving particle peculiar 
velocity, and x is the comoving position. Throughout, 
dots refer to derivatives with respect to proper time t. 
Perturbations in the scalar field give rise to the fifth force 
on the right-hand side in the equation above. 
The comoving potential satisfies the Poisson equation: 
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with p as the comoving DM particle density. Here and 
throughout, an overline indicates a spatial average. 

At the present epoch, the field behaves as a cosmolog- 
ical constant, so the potential term in Eq. (TT] dominates 
and has a value 



V((f>b,o) = ^AOpcrit- 



(4) 



At early times, the coupling to matter dominates the 
scalar field equation of motion, and Eq. (JTJ) reduces to 
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which we use to set the initial condition for fa. 
The DM particle also has a field-dependent mass 
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The Friedmann equation, neglecting radiation, curva- 
ture, and baryonic terms, becomes 



— — H n fl r 



5 'mO 

a I 0b,o 



b -a- 3 



8ttG 



Kfi — 



(7) 

To simulate a comparative ACDM cosmology, we fixed 
(/>b to the value in Eq. ^ and prevented any interactions 
between the field and particles. 

1.2. Numerical Methods 
We study a general power-law potential: 

V(<f>)=K/<f> a , (8) 

where we are free to choose the constants K and a. Based 
on the comments made by Farrar and Peebles and our 
analysis in SR08, we chose two combinations of param- 
eters. These were selected for behaviors that were sig- 
nificantly different from standard cosmologies, but not 
drastic enough to rule them out with current observa- 
tional contraints. Table [1] lists the parameter values, the 
guessed initial field value at our simulation initial red- 
shift, and the field value today as calculated from Eq. ([1]). 
For our simulations we chose FLASH v2.5, an adaptive- 
mesh refi nement (AMR) cod e for astrophysics and cos- 
mology (|Frvxell et all l2000fl . FLASH solves the N- 
body potentia l problem wi th a particle-mesh multigrid 
FFT method fjRicker|[20M) . FLASH uses cloud-in-cell 
mappi ng for interpolating between the mesh and par- 
ticles (Hockney & East wood! 119881) and a second-order 



leapfrog integration scheme for variable-timestep particle 
advancement. We modified the standard FLASH code by 
adding the additional drag and force terms from Eq. $2$. 
At each timestep, the particle mass is updated according 
to Eq. ([6]). We calculate the scale factor and scalar field 
value in-code by numerically solving Eqs. ((7|) and ([lj, 
respectively. For a more detailed explanation of solving 
the scalar field equation , see SR08. 

For all calculations, we used concordance parame- 
ter values of fi m o — 0.26, fl\o — 0.74, and Hq = 
lOO/i = 71 km s~ Mpc" 1 . All runs took place in a 
three-dimensional box measuring 128 hr 1 Mpc per side 
with 256 3 particles. For each model, including a ACDM 
reference, we performed 10 simulations with 512 zones 
per side and an additional 4 simulations with 1024 zones 
per side to study lower mass ranges. There was no re- 
finement of grid spacing. All simulations used the same 
initial conditions: unperturbed particle positions were 
situated on a grid, and the initial velocities and positions 
were perturbed using Guassian fluctuations normalized 
to erg = 0.751. We assumed P(k) fro m a ACDM cosmol- 
ogy. We used the GRAFIC2 code (jBertschingerll2001h 
to generate these initial conditions. All computations 
started at a redshift of z = 56.8. 

We used a friends-of-friends (FOF) routine to find ha- 
los. This algorithm builds lists of all particles that are 
within a certain maximum distance of their neighbors. 
For all results, we chose a linking distance of 1/5 of the 
unperturbed particle spacing, which is 500 h~ x kpc. At 
every analysis redshift, we calculated the minimum re- 
so lvable halo particle count according to the prescription 
in lLukicet atl (|2007h : 
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where n g and n p are the number of zones and particles 
per side, respectively. We chose an overdensity factor of 
A = 200. To examine the halo mass function, we cor- 
re cted the halo FOF p article counts by the factor given 
in lWarren et alJ (|2006h : 
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2. THE HALO MASS FUNCTION 

Figure [T] shows the relative mass function from redshift 
z = to z = 1.5 for the two models listed above. Both 
of these are compared to the ACDM simulation mass 
function. We will use the frequency density definition of 
the mass function, 



F(M, z) = log 
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dlogM' 
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so that for model i we may define the relative mass func- 
tion as 



RF(M, z)i = log 
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log 
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ACDM 



We analyzed relative mass functions to reduce any sys- 
tematic errors in the simulations, including those due to 
small box effects, such as missing tidal forces. We binned 
our distributions into 10 fixed logarithmic intervals from 
10 115 to 10 15 Ii~ 1 Mq. We only display values in bins 
for which we have complete data (i.e. the bin does not 



contain the minimum resolvable mass). The uncertain- 
ties shown are obtained by summing in quadrature the 
individual statistical counting errors in the interacting 
and ACDM cases. We see that at the present epoch, 
both models produce a greater number of the most mas- 
sive halos while underproducing low-mass objects. At 
higher redshifts, both models produce greater numbers 
of all objects. Since the relative mass function does not 
remain constant with redshift, it is distinguishable from 
a concordance cosmology with different fundamental pa- 
rameters. 

To determine the significance of these mass function 
differences, we performed a \ 2 likelihood test at each 
redshift. For two indepedent frequency distributions R 
and S 1 , 

-r-^ (Si — Ri) 



X 
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(13) 



where the sum takes place over all bins and the num- 
ber of degrees of freedom is the total number of non-zero 
bins. Figure [5] shows the probability at each redshift that 
the frequency distributions from the interacting cases are 
consistent with the ACDM case. We see that for al- 
most all redshifts, the probability is exceedingly small, 
indicating that these coupled DM-DE models are well 
distinguished from ACDM cosmology with this sample 
of objects. As a function of redshift, the probability for 
both models generally decreases. For model A at z = 0.5, 
however, these distributions are not well separated, since 
this is the redshift at which the mass functions come 
closest. Model B remains indistinct from ACDM until 
z = 0.5, at which point the separation becomes progres- 
sively more evident. The probabilities for both models 
drop below 10~ 10 at redshifts z = 0.75 and z = 1.5 and 
are not shown at higher redshifts. 

To determine the feasibility of using a survey 
like DES to detect these models, we need to know 
the expected frequency distribution of observed 
objects in the survey, since the fundamental uncer- 
tainty arises from counting errors. While the exact 
selection function for DES has not yet been deter- 
mined, the survey is expected to capture ~ 10,000 
obj ects of mass > 10 13 5 h~ 1 M P] out to redshift z ~ 
1.5 dThe Dark Energy Survey. White Paper submitted to 
2005). Our combined simulations produced roughly 
this many objects in this mass and redshift range. 
Figure [2] shows the \ 2 probability when only considering 
mass bins above 10 135 Ii~ 1 Mq. With this cut the 
probability suffers; however, we maintain the general 
trend of increased disparity with higher redshift. For 
both models, the greatest deviations occur at z = 0.75 
and z = 1.5. The probabilities at these redshifts are 
below the common significance threshold and indicate 
that DES is capable of detecting these models. At lower 
redshifts, there is not a significant difference in the 
number of high-mass objects. Between these redshifts, 
the uncertainties in the largest mass bins are too large to 
make a confident distinction between the distributions. 

It is possible to observe similar deviations in the mass 
function by changing the fundamental parameters of a 
ACDM cosmology, such as ft m o an d cr§. We could not 
find any combination of fundamental parameters that re- 
produced these relative mass functions for all mass bins 
and redshifts. However, when examining masses within 



the DES limit, there are degeneracies. For example, 
Figure [T] shows the re l ative mass function, obtained us- 
ing the IWarren et al.l (|2006f ) analytic mass function, of 
two ACDM cosmologies with different values of cr§. In 
this case, we compared a og of 0.751, which we used in 
our simulations, to a cosmology with erg = 0.775. We 
chose this value to mimic the entire relative mass func- 
tion at high redshift, but it does not capture the counts 
of lower-mass objects at low redshift. However, above 
10 135 h~ 1 M & , the uncertainties are large enough to per- 
mit this modified erg to adequately fit the data. Thus, 
DES alone may not be able to distinguish between a uni- 
verse with interacting DM-DE and a universe with higher 

While these two models are easily distinguishable from 
ACDM cosmology, they are more difficult to differenti- 
ate from each other. Figure [3] shows the % 2 probability 
when comparing these two models, both for all masses 
and DES-accessible masses. When including all masses, 
the probabilities are significant at redshifts z = 0.0, 1.0, 
and 1.5. However, when examining only masses available 
to DES, the probabilities maintain a roughly constant 
performance and never reach a significant level. 

3. CONCLUSIONS 

We have found that coupling dark matter to a dark 
energy scalar field produces significantly different mass 
functions at redshifts as high as z = 1.5 relative to a 
ACDM cosmology with the same set of fundamental pa- 
rameters. This difference in the mass function follows 
from our analysis in SR08: an additional fifth force and 
a reduced particle Hubble drag lead to more structures 
than in ACDM cosmologies at early times, and at late 
times will cause an overabundance of high-mass objects 
and a subsquent reduction in low-mass cluster counts. 
By examining the mass function, we have developed a 
simple way of distinguishing these models. This analysis 
allows us to discover ways of further constraining differ- 
ent parameters of DM-DE coupling. 

We have found that the statistical uncertainties in the 
mass function do not prevent the Dark Energy Survey 
from detecting this form of coupled dark matter and dark 
energ y. Once the selection function for DES is known, a 



the spJa^Efbcm^iBaBlmFxHccan be applied and a more de- 
tailed study will need to take place. However, we have 
found that DES alone will have difficulty differentiating 
among different sets of parameters that control the cou- 
pling. 

At high masses the statistical uncertainties may pre- 
vent DES from distinguishing between coupled DM-DE 
and ACDM cosmologies with different values of erg ■ We 
can overcome this dege neracy in seve ral ways. First, mis- 
sions such as Planck (jTauberl [ 20 00) can independently 
constrain og and £l m o- If DES prefers a higher value 
of erg through the mass function, this may be explained 
by interacting dark matter and dark energy. Secondly, 
more detailed measurements of the halo m ass function 
from projects such as LSST (|Sweenevll2006| ) will tightly 
constrain the mass function at multiple redshifts. Also, 
DES itself may detect more clusters than our estimated 
10,000. A universe with coupled DM-DE will then pro- 
duce an apparently redshift-dependent erg. 

It will also be necessary to compare these mass func- 
tions to those produced by modified General Relativity 
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Fig. 1. — Relative mass functions for the power-law potentials.. The plots are labeled as in Table [T] Error bars are at one standard 
deviation and come from statistical counting uncertainties. The solid curve is the relative mass function of two ACDM cosmologies with 
ct 8 = 0.775 and cr 8 = 0.751. 

(such as those found in lStabenau fc Ja in (2006)), as both 
theories modify the Poisson equation, and hence can in 
principle have similar effects. 
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Fig. 2. — Probability for the chi-square statistic as a function 
of redshift for model A (top) and model B (bottom) compared to 
ACDM. The solid lines are from including all resolvable masses, and 
the dashed lines are from only considering objects with MpoF > 
10 13.5 k -l M 
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Fig. 3. — Probability for the chi-square statistic as a function of 
redshift for model A compared to model B. The solid line is from 
including all resolvable masses, and the dashed line is from only 
considering objects with MpQF > 10 13 ' 5 h~ 1 MQ. 
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